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Abstract. We present a 23 ks XMM-Newton observation of the Ultraluminous Infrared Galaxy (ULIRG) Mrk 273. The hard 
X-ray spectrum can be modeled by a highly absorbed (~ 7 x 10^^ cm"^) power law plus an Fe Ka emission line. The iron 
line (detected at more than 99% c.l.) is broad (cr = 0.26^q keV), suggesting possible superposition of a neutral iron line 
at 6.4 keV, and a blend of ionized iron lines from Fe XXV and Fe XXVI. Given the relatively short exposure, the three line 
components can not be singularly resolved with high statistical significance: the neutral component is detected at ~ 2.5cr and the 
Fe XXV line at ~ 2it c.l., while for the Fe XXVI line we can only estimate an upper limit. The broad band spectrum requires, 
in addition to a highly absorbed power law, at least three collisionally ionized plasma components, which may be associated 
with star-forming regions. The temperatures of the three plasmas are about 0.3, 0.8 and 6 keV, where the highest of the three is 
sufficient to produce ionized iron emission lines. An alternative interpretation for the origin of the soft emission might also be 
given in terms of reflection off some photoionized gas, as has been observed in a number of nearby Compton-thick Seyfert 2 
galaxies (e.g. NGC 1068, Circinus, Mrk 3, NGC 4945). A hot gas, photoionized by the primary, continuum can also produce 
ionized iron lines. Unfortunately, given the limited statistics and the lack of high resolution spectroscopy, it is not possible 
to distinguish between the two models investigated. We further compare the XMM-A'evvton findings with the Chandra data 
obtaining consistent spectral results. The absorption corrected hard X-ray luminosity of Mrk 273 is Li-iokev ~ 7 x lO'*^ erg s"', 
corresponding to ~ 0.2% of the far-IR luminosity, similar to typical values found in pure starbursts. The thermal contribution to 
the soft X-ray luminosity is approximately 0.2 - 0.7 x 10*^ erg s"' , comparable to those found in NGC 6240 and other starburst 
dominated ULIRGs. We also analyze the XMM-Newton spectrum of Mrk 273x, an unabsorbed Seyfert 2 galaxy at redshift 
z = 0.458, which lies in the field of view of Mrk 273. 

Key words, galaxies: individual: Mrk 273 - galaxies: individual: Mrk 273x - galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

The InfraRed Astronomical Satellite (IRAS) has detected a 
very large number of galaxies in the local universe (z < 0.3) 
which exhibit extraordinarily high infrared luminosity (Lir > 
10^^ Lo for Ho = 75kms 'Mpc ^ see Sanders & Mirabel 
1996 for a review). The number density of these so called 
Ultraluminous InfraRed Galaxies (ULIRGs) exceeds that of 
optically selected Seyfert ga laxies and QSOs v vith compara- 



id q: . 

ble bolometric luminosities ('Soifer e t alJll987t ISanders et alJ 
[l988ab; Sanders & Mirabel 1996) by a factor of ~ 1.5-2 
( ISanders et al J 19991) . The bulk of luminosity in these sources is 
infrared (IR) emission from warm dust. Only two dust-heating 
mechanisms are capable of producing such an extraordinary IR 
luminosity: one involves the presence of a strong starburst re- 
gion, and the other a dust-enshrouded AGN. However it is still 
not clear what the relative contribution of each component to 
their bolometric luminosity should be. 



Spectroscopic surveys of samples of ULIRGs carried out 
with the Infrared Space Observatory (ISO) have revealed that 
about 80% of such objects are predominantly powered by star 
formation, but the fraction of AGN-powered objects increases 
with luminosity cGenzel et al. 1998: Lutz et al. 1998) . 

X-ray observations are a fundamental tool to probe the 
highly obscured innermost regions of these objects, and there- 
fore potentially unveil the physical processes at work. In fact 
a considerable portion of ULIRGs has been found to con- 



tain a hard X-rav source, highlv absorbed bv a molecular 
torus, which indicates the presence of a hidden AGN liMitsudal 


1995; 


Brandt et al. 


1997; Kii et al. 


119971: IVienati etalJll999l: 


Ptak et alJ2003l;lBraito et alJ2003l 


Franceschini et alJ2003t). 



Send offprint requests to: Italo Balestra 
e-mail: balestraSmpe . mpg . de 



Mrk 273 is a well studied ULIRG at a redshift z = 0.03778 
with a Seyfert 2 nucleus ( Koski 1978: Sanders et al. 1988al). 
It pre sents evidences for strong star formation ( iGoldader et alJ 
11995V such as strong Pol ycylic Aro matic Hydrocarbons (PAH) 
features in the near-IR (iGenzel et a l. 1998) and extended (~ 
370 pc) radio emission, punctuated by a number of com- 
pact sources identifiable as supernovae or supernova rem- 
nants (iCarilli & Tavlor. .200(1 .Bondi et aL .20051) . Like many 
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ULIRGs, Mrk 273 presents clear indication of an ongoing 
merging process such as a long tidal tail and a double riucleu s 
llKnanen et alJl997HSoifer et alJ200(ilCarim & Tavloil20 00Y 

ASCA X-ray observations of Mrk 273 showed evidence 
for a highly absorbed (~ 4 x 10^^ cm"^) hard component above 
about 3 keV together with soft thermal emission llTurner et alJ 
ll 9971 ll 99Rt iTwasawalfl 999.) . These findings were indeed con- 
firmed by Chandra, which revealed a compact hard X-ray nu- 
cleus inside a much more extended soft halo llXiaetalJl20()i 
hereafter Paper I). 

In this paper we will discuss the analysis of the XMM- 
Newton data of Mrk 273 and compare them with the previous 
Chandra observation (Paper I). The data analysis and details 
of the 'KMM.-Newton observation are described in Section 2, 
together with our updated reduction of the Chandra data. In 
Section 3 the spectral fitting results are presented. A com- 
parison with the Chandra observation is given in Section 4. 
Our results are then discussed in Section 5 and a summary of 
the XMM-Newton results on Mrk 273 is given in Section 6. 
In Appendix A we discuss the XMM-Afewfon spectrum of 
Mrk 273x, an unabsorbed Seyfert 2 observed serendipitously 
in the field of view of Mrk 273. 

Throughout this paper we use the following cosmological 
parameters: //q = 70 km s"' Mpc"', Aq - 0.7 and q'o = 0. 
Within the adopted cosmology 1" coiTesponds to 0.75 kpc at 
theredshiftofMrk273. 

2. Observations and Data Reduction 

2.1. XMM-Newton 

Mrk 273 was observed with XMM-Newton jjansen et al.l200ll) 
for 23 ks on 2002 May 07 (revolution 0441). During this 
time the EPIC-PN (Struder et al. 200D and MOS (MOSl and 
MOS2; Turner et al. 2001) cameras, as well as the Optical 
Monitor (OM; Mason et al. 2001) and the Reflection Grating 
Spectrometers (RGSl and RGS2; den Herder et al. 2001) col- 
lected data. The EPIC-PN and MOS cameras were operated in 
full-frame mode and utilized the thick filter. 

The Observation Data Files were processed to produce cal- 
ibrated event lists using the XMM-A^ewfo« Science Analysis 
System (SAS v6.1.0). Unwanted hot, dead, or flickering 
pixels were removed as were events due to electronic noise. 
Event energies were corrected for charge-transfer losses, and 
EPIC response matrices were generated using the SAS tasks 
ARFGEN and RMFGEN. Light curves were extracted from these 
event lists to search for periods of high background flaring. 
Background flaring was negligible. The total good exposure 
times selected for the PN and MOS were 18 and 22 ks, re- 
spectively. 

The source plus background photons were extracted from 
a circular region with a radius of 50", and the background was 
selected from an off-source region with a radius of 70" and 
appropriately scaled to the source region. Single and double 
events were selected for the PN detector, and single-quadruple 
events were selected for the MOS. The resulting PHA files 
were grouped with a minimum of 20 counts per bin. Pile-up 
effects were determined to be negligible. 



The RGS were operated in standard SpectroH-Q mode. The 
first-order RGS spectra were extracted using the SAS task 
RGSPROC, and the response matrices were generated using 
RGSRMFGEN. Unfortunately, due to the low signal-to-noise the 
RGS data were not fruitful. 

The OM was operated in imaging mode for the entire obser- 
vation. Fourteen images were taken in three filters: 5 in UVW\ 
(245 - 320 nm), 9 in UVM2 (205 - 245 nm), and 1 in UVW2 
(180 - 225 nm). The average exposure time (800 s) was short 
for such a faint source and magnitudes could not be calculated 
for all UVM2 images. The average apparent magnitude in each 
filter was UVWl = 16.40 + 0.07, UVM2 = 16.51 + 0.08, and 
UVW2 - 16. 14+0.07, somewhat larger than the ground-based 
measured U (309 - 373 nm) magnitude ni jj — 15.1 + 0.2 for 
the whole galaxy ( Surace & S anders 2000), but not inconsis- 
tent when accounting for the different aperture sizes used. 

2.2. Chandra 

Chandra observed Mrk 273 on April 2000 with the Advanced 
CCD Imaging Spectrometer (ACIS-S3 Back Illuminated chip) 
for about 47 ks. The data relevant to this observation were 
already published in Paper I, where they were processed fol- 
lowing the Standard Data Processing available at that time 
(R4CU5UPD13.2, January 2001). 

Here, data are reprocessed using much more recent versions 
of the Chandra Interactive Analysis of Observations software 
(CIAO 3.2) and Chandra Calibration Database (CALDB 3.0.0). 
The most effective improvements introduced concern the more 
accurate background subtraction and more careful computation 
of the effective areas and response matrices. 

We start processing data from the level=l event file. We 
apply the recently released, time-dependent gain correction \ 
which is necessary to adjust the "effective gains", which have 
been drifting with time due to increasing charge transfer in- 
efficiency. Since the observation was taken in the VFAINT 
mode, we run the tool acis_process_events to flag prob- 
able background events using all the information of the pulse 
heights in a 5 x 5 event island to help distinguishing between 
good X-ray events and bad events that are most likely associ- 
ated with cosmic rays. With this procedure, the ACIS particle 
background can be reduced significantly compared to the stan- 
dard grade selection^. Real X-ray photons are practically not 
affected by such cleaning (only about 2% of them are rejected, 
independently of the energy band, provided there is no pileup). 

The data are filtered to include only the standard event 
grades 0, 2, 3, 4 and 6. We finally filter time intervals with high 
background by performing a 3cr clipping of the background 
level using the script analyze_ltcrv-'. This yield an effec- 
tive exposure time of about 43 ks for the AC1S-S3 chip in the 
energy range 0.3 - 10 keV. 

As in Paper I, we separate the spectral analysis of the nu- 
clear emission from the extended soft X-ray halo. The spec- 
trum of the nuclear region of the source is extracted from a 

' http ://asc . harvard. edu/ ciao/threads/acistimegain/ 

^ http://asc.harvard.edu/cal/Links/Acis/acis/CaLprods/vfbkgmd/ 

^ http://cxc.harvard.edu /ciao/tlireads/filterjtcrv/ 
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Fig. 1. Adaptively smoothed 0.3 - 10 keV Chandra image of 
Mrk 273 with a smoothing scale of 2 pixels (w 1"). The circle 
has a radius of 10" and it is centered at the peak of the hard X- 
ray emission. The ellipse has a size of 42" x 25" and a position 
angle of 13°. The bright source visible in the upper left corner 
is Mrk 273x. 



Table 1. XMM-Newton and Chandra net number of counts 
(background subtracted) collected by each instrument in the 
specified extraction region and spectral range. 



Instrument 


Extraction 
Region 


Energy 
[keV] 


Net Number 
of Counts 


XMM-Newton 










PN 


50" 


0.3- 


- 10 


1680 


MOS 1 


50" 


0.3- 


- 10 


590 


MOS 2 


50" 


0.3- 


- 10 


520 


Chandra 










ACIS-S 


10" 


0.3- 


- 10 


1990 


ACIS-S 


ellipse 


0.3- 


-2 


530 



Throughout our spectral an alysis we fix the abundanc es of 
the elements to the solar values dAnders & Grevessell989l) and 
we model Galactic absorption with the photoelectric absorption 
model PHABS in XSPEC, where we fixed the absorbing col- 
umn density to Ng = 1.09 x lO^^ cm ^ (HEASARC W3nH4). 
Line energies are given in the rest frame of the source. 



circular region of radius 10" centered at the peak of the hard 
X-ray emission, while, for the extended soft halo, we choose 
the region lying between the inner circle and the ellipse shown 
in Figure [0 in order to have the best possible signal-to-noise 
ratio. 

We use the events included in each of the extraction re- 
gions defined above to produce a spectrum (PHA) file, which 
we grouped with a minimum of 20 counts per bin, unless other- 
wise specified. The background is obtained from a source free 
circular region of radius 50" located on the same chip. The 
background file is then scaled to the source file by the ratio 
of the geometrical area. The response matrices and the ancil- 
lary response matrices of each spectrum are computed with 
acisspec for the same regions from which the spectra are 
extracted. 

The X-ray spectral analysis has been performed using 
XsPEC version 11.3.1. 

The reported errors on model parameters are at the 90% 
confidence level for one interesting parameter (Ax^ - 2.71). 

3. Spectral Analysis 

In Table we give a comparison between the XMM-Newton 
and Chandra data sets. The Chandra ACIS-S spectrum of the 
nuclear region (inner 10") of Mrk 273 is dominated by the 
background above 8 keV, while the spectrum of the extended 
soft X-ray halo is source dominated only between 0.5 and 0.9 
keV. Both the EPIC-PN and the two MOS spectra are source 
dominated between 0.3 and 10 keV, therefore the entire energy 
range can be analyzed with relatively good signal-to-noise. 
XMM-Newton therefore provides some further constraints on 
the shape of the continuum between 8 and 10 keV. 

The extracted PN and ACIS-S light curves do not show any 
significant (within 30%) flux or spectral variability. 



3. 1. The complex Fe Ka line 

In order to study the properties of the Fe Ko' line we concentrate 
our spectral analysis only on the hard X-ray (4-10 keV) spec- 
trum of the EPIC-PN. In this energy range the continuum can 
not be reproduced by a simple unabsorbed power law, since this 
model gives a negative spectral index and a poor fit. Allowing 
for extra absorption at the redshift of the source (ZPHABS in 
XsPEC), we find a much better fit result (xj ^ 1.2), but, given 
the limited statistics available, if both F and Nh are left free to 
vary, we find very large uncertainties (e.g. F - 2.4;^^^). Fixing 
the spectral index to the "canonical" value F = 1.9 allows us 
to measure a column density Nh = (6.2 + 1.5) x 10^^ cm"^ in- 
trinsic to the source. As displayed in Figure|2] the fit is still not 
completely satisfactory (x^/dof = 24.5/20), mainly because 
of the presence of a feature in the spectrum at about 6.4 keV, 
which can be identified with an Fe Ka emission line. A sig- 
nificant improvement in the fit is obtained when a single broad 
Gaussian is added to the model (x^/dof - 15.2/17). The line 
centroid is at 6.51 + 0.19 keV and the line is found to be broad 
(cr = 0.26+J}^^ keV), with an equivalent width EW = 560+^^0 
eV. 

Such a broad iron line (FWHM- 29 000 km s ') could, in 
principle, originate from the accretion disc, since the primary 
radiation is obscured only up to about 4-5 keV and therefore 
should only partly affect the photons at the iron line energy. We 
try to model the iron line profile with the DISKLINE model 
in XsPEC for a Schwarzschild black hole, which parameterizes 
the radial emissivity of the disc as a power law (i.e. oc r"?). We 
fix the line rest energy to 6.4 keV, the emissivity of the disc q to 
-2, the inner and outer radii respectively to 6 and to 1000 rg and 
the inclination angle to 30°. We find an equivalent width EW = 
520 + 300 eV and a fit result statistically equivalent to the one 

* http;//heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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Energy (keV) 

Fig. 2. The 4.0- 10 keV EPIC-PN spectrum fitted with a simple 
absorbed power law model (F is fixed to 1.9). Residuals are 
found at Fe Ka emission line energy (at about 6.5 keV). 



obtained with a sing le broad Gaussian (x^/dof = 17.1/19). 
Unfortunately, the statistics are too poor to derive any physical 
property of the disc. 

The hypothesis of a broad iron line originating from the 
disc can not be ruled out, but the high column density, which is 
obscuring the AGN up to about 4 keV, certainly makes the de- 
tection of the redshifted wing of a relativistic iron line profile 
much less feasible. A more likely explanation for the promi- 
nent iron line feature observed may be given in terms of a 
complex profile due to the superposition of, at least, two un- 
resolved line components: a neutral iron line at 6.4 keV, prob- 
ably originating from a Compton-thick torus, and an ionized 
(Fe XXV) iron line at 6.7 keV. Indeed, the latter could be pro- 
duced either by some diffuse coUisionally ionized plasma at 
temperature kT > 5 keV, as observed in at least another two 
ULIRGs (namely NGC 6240, Boiler et al. 2003; Netzer et al. 
2005 and Arp 220, Iwasawa et al. 2005), or by reflection off 
some photoionized gas surrounding the nucleus, as often seen 
in Compton-thick Seyfert galaxies (e.g. Circinus, Sambruna 
et al. 2001; NGC 1068, Kinkhabwala et al. 2002; NGC 4945, 
Done et al. 2003; Mrk 3, Bianchi et al. 2005). 

When the line profile is modeled with two narrow (cr - Q 
eV) Gaussian lines with energies fixed at 6.4 and 6.7 keV the 
result is statistically equivalent to a single broad line model 
(X^/dof = 15.5/18). In this case we measure an equivalent 
width EW = 165![^|| eV and £W = 120 + 95 eV for the neutral 
and for the ionized component respectively. Finally, we find no 
improvement in the fit for the addition of a third ionized iron 
line at 6.97 keV from Fe XXVI. However upper limits on the 
flux (< 0.18 X 10 erg cm"^ s ') and the equivalent width 
(< 85 eV) of this line can be derived. 

In the following we discuss two distinct models capable of 
reproducing the broad band spectrum as well as the possible 
ionized iron lines: one through thermal emission, therefore in- 
volving the presence of a hot gas associated with a starburst, 
and the other one through reflection of the primary continuum 
by some optically thin photoionized gas. 



Table 2. XMM-Newton EPlC-PN spectral fitting results when 
the thermal emission model is applied (2 APEC: two temper- 
ature model; 3 APEC: three temperature model). The last col- 
umn displays the result of the combined PNh-MOS fit with the 
3 APEC model. Fluxes are given in units of 10 erg cm"^ 
s luminosities in units of 10'*^ erg s"'; both are corrected for 
absorption. * denotes a fixed parameter. 





2 APEC 


3 APEC 


3 APEC 




PN 


PN 


PN + MOS 


Apec(l) 








kT [keV] 


^•'-1.1 


"■^^-2.0 




/o.3-10keV 


1.40 ±0.15 


1.36 ±0.18 


1.43 ±0.14 


Apec(2) 








kT [keV] 


"•""-0.04 


0.75 ± 0.07 


0.78 ± 0.08 


/ 0.3-10 keV 


^ ^ -^4-0 1 1 

63 

"■"■^-0.07 




'-'■■^-^-0.15 


47 

"■^'-0.15 


Apec(3) 

kl [keVJ 

/ 0.3-10 kcV 




n 08 

23+° 

"•^-'-0.10 


n T 1 +0.09 
"•-'^-0.07 

23+" '^ 

"•^-^-0.09 


Power-law 










68!]^ 


59+16 


59+15 

"^-14 


r 


1.9* 


1.9* 


1.9* 


/o.3-10keV 


44+21 

^^-13 


44+27 


46!!? 


FeKa 








EW [eV] 

/o.3-10kcV 


0.48!«;« 


91 Q + 198 

49+0 

"•^^-0.33 


130 ± 120 
0.31 ±0.29 


Total 








/o.3-10keV 


46.6 


46.6 


48.7 


-i^0.3-10kcV 


1.54 


1.54 


1.66 




84.2/80 


63.8/78 


130.0/134 



3.2. The thermal emission model 

At first we concentrate our analysis on the EPIC-PN spectrum. 
The better statistics, especially at high energies, permit more 
robust constraints of the model parameters. Once the best fit 
for the PN is found, we perform a combined fit of the PN and 
the two MOS, applying the same model, in order to further 
improve the statistics. 

When the 0.3 - 10 keV EPIC-PN spectrum is modeled wifli 
a simple absorbed (Galactic plus intrinsic) power law, we ob- 
tain a statistically unacceptable fit (x^ > 4). Strong residuals 
remain primarily above 4 keV, but also between 0.6 and 1 keV 
(Fig.EJ. 

As a starting point for our spectral modelling we adopt the 
best fit absorbed power law plus Gaussian model for the 4-10 
keV band (see Section 3.1) and extrapolate to 0.3 keV. A very 
prominent soft excess is clearly noticeable. Since we intend to 
investigate both the origin of the soft excess and of the possi- 
ble ionized iron, in the following we fix the parameters of the 
Gaussian profile to E = 6.4 keV and cr - eV. This will be 
considered our baseline model. 

Diffuse, optically thin, thermal plasma emission is often 
observed in Starburst galaxies and ULIRGs. It is typically as- 
sociated with star-forming regions and shock-heated gas in 
merging systems. We try to model the soft, extended. X-ray 
emission with a thermal spectrum from a coUisionally ion- 
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Fig. 3. A simple absorbed power law model applied to the 
EPIC-PN spectrum of Mrk 273. Strong residuals (data/model) 
are clearly present above 4 keV and also between 0.6 and 1 
keV. 

ized plasma, using the APEC model dSmith et alJl200ll) in 
XSPEC. A single temperature plasma component gives a con- 
siderable improvement to the fit with respect to the baseline 
model (x^/dof = 259/82), but can not simultaneously account 
for the observed peak at about 0.7 - 0.9 keV and the possible 
ionized iron lines. 

Therefore, we add to the baseline model two collisionally 
ionized plasma components of different temperatures. The free 
parameters are now the temperature and the normalization of 
each plasma component, the column density of the absorber at 
the redshift of the source, the normalization of the power law 
and the normalization of the Gaussian line. The best fit param- 
eters for this model are given in Table|2] A plasma temperature 
of ~ 0.7 keV is needed in order to model the peak in the emis- 
sion at low energies, while a much hotter (~ 5 keV) plasma 
is required to fit part of the high energy continuum and also 
to give rise to the ionized iron line at 6.7 keV. The second, 
hotter, plasma is required with a high statistical significance 
(Ax^ ~ 175 for the addition of 2 free parameters). However the 
fit is still not completely satisfactory (x^/dof - 84.2/80) due 
to the presence of some excess below about 0.6 keV. 

The addition of a third, lower temperature, plasma compo- 
nent provides a further improvement to the fit (A^^ - 20 for the 
addition of 2 free parameters). In this case we find an undoubt- 
edly acceptable result (x^/dof - 63.8/78). The lowest tem- 
perature component is required with a probability of more than 
99.99% according to an F-test. Its temperature (kT - 0.26+qo^ 
keV) is consistent with the one derived from the spectrum of 
the extended hot gas halo (kT = 0.37 + 0.04 keV) observed in 
the Chandra image (see Section 4.2). 

The best fit to the 0.3 - 10 keV PN spectrum is therefore 
obtained using three thermal plasma components with differ- 
ent temperatures, a highly absorbed power law and a neutral Fe 
Ka line (Figure 1^. According to this model the emission be- 
low about 4 keV would be purely thermal. We will discuss the 
possible implications of these new findings in Section 5. 

Simultaneously fitting the PN and the two MOS data results 
in an equivalent fit, as shown in Figure^] All the parameters are 
consistent with those obtained using the PN only (see Table |3. 
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Fig. 4. (Upper Panel) The EPIC-PN spectrum of Mrk 273 fit- 
ted with three different temperature APEC models, a highly 
absorbed power law and a neutral Fe Ka line. Each model com- 
ponent is displayed. (Lower Panel) The EPIC-PN and the two 
MOS combined fit with the same model. The closed and open 
circles correspond to the two MOS data sets. The best fit pa- 
rameters are given in Table |2] 



It is worth noting that, combining the PN and the MOS, we find 
a temperature kT = 0.31 ^{5°^ keV for the coolest thermal com- 
ponent, which is in very good agreement with the one measured 
by Chandra for the hot extended halo. 

Finally, we remark that there is no need for further absorp- 
tion at the redshift of the source for any of the three thermal 
components. When left free to vary the absorbing column den- 
sity of each component is consistent with the Galactic value 
and do not improve the fit (Ax^ < 1 for the addition of one free 
parameter). 

3.3. The photoionized gas model 

An alternative interpretation for the origin of the soft X-ray 
continuum and the ionized iron line observed in the XMM- 
Newton spectrum is scattering of the primary radiation by some 
photoionized, Compton-thin, hot gas surrounding the AGN. 
Ionized Fe Ka emission lines from H- and He-like iron ions 
are expected to be produced through recombination and reso- 
nant scattering jfiianchi & Mati2002h . 

The analysis of the Chandra image in Paper I revealed that 
the emission below 1 keV is extended on a region of radius 
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~ 7.5 kpc surrounding the nucleus and, partly, on a much larger 
region (~ 45 kpc) encompassing the tidal tail. A pure reflec- 
tion model can not account for such extended X-ray emission, 
which must be of thermal origin. If we model the 0.3 - 10 keV 
PN continuum with a highly absorbed (primary) plus a less ab- 
sorbed (reflected) power law having the same spectral index, 
we find an unacceptable fit (x^ ^ 2.6). Therefore we added to 
the model two narrow Gaussian lines, at 6.4 keV from neutral 
iron and at 6.7 keV from Fe XXV, plus a hot coUisionally ion- 
ized plasma component. In this case the fit result is acceptable 
(X^/dof - 66.9/79). However, there is no need for an addi- 
tional lower temperature thermal plasma component, which is 
clearly detected by Chandra. 

This apparent discrepancy between the two data sets may 
be reconciled simply by taking into account the absorption of 
the secondary power law, which is apparently negligible in the 
XMM-A^ewto« specti'um, but ~ 7 x 10^" cm"^ in the Chandra 
spectrum of the "nuclear" region (see Section 4.1). Therefore 
we fix the column density of the material obscuring the re- 
flected power law to the value measured by Chandra. The free 
parameters of the model are: the temperature and the normal- 
ization of the thermal component; the column density of the 
absorber at the redshift of the source, which is obscuring only 
the primary power law; the normalizations of the two power 
law components; and the normalizations of the two Gaussian 
lines. 

When this model is applied to the PN spectrum the fit result 
is still not completely satisfactory (jc^/dof = 84.9/80), mostly 
due to some excess still present below 0.6 keV. The best fit 
parameters for this model are given in Table|3] 

The addition of a second, lower temperature, plasma com- 
ponent is a considerable improvement to the fit (A^^ ^ 20 for 
the addition of 2 free parameteis, / do f - 64.0/78). Figure|5l 
shows the best fit obtained for this model. The cooler ther- 
mal plasma component is required with a probability of more 
than 99.99% according to an F-test. The temperature found 
kT - 0.26;^Q Q3 keV is consistent, within the uncertainties, with 
that measured by Chandra for the extended hot gas halo. 

The flux of the reflected power law amounts to 2 - 6% of 
the primary component. From this ratio we can estimate the 
column density of the reflecting material. Assuming a cover- 
ing factor of 0.5, the column density of the photoionized gas 
would be approximately 10^^ cm"^. The equivalent widths of 
the ionized iron lines expected to be produced by a hot gas with 
such a density are consistent with the measured values within 
the uncertainties. 

The best fit model obtained here is similar to the one pre- 
sented in Paper I, with the only noticeable differences being a 
slightly larger absorbing column density obscuring the primary 
radiation and a slightly smaUer temperature of the cooler ther- 
mal component. 

4. Comparison with Chandra 

For a detailed spatial analysis of the Chandra X-ray image of 
Mrk 273 we refer the reader to Paper I. Here, as in Paper I, we 
extract separately the spectra of the nuclear region and that of 



Table 3. XMM-Newton EPIC-PN spectral fitting results using 
the photoionized gas model described in the text. Fluxes are 
given in units of lO"'-' erg cm~^ s"\ luminosities in units of 
10"*^ erg s"' ; both are corrected for absorption. * denotes a fixed 
parameter 
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Fig. 5. The EPIC-PN spectrum fitted with two thermal plasma 
components, a highly absorbed and a less absorbed power law 
with the same spectral index, plus a narrow Fe Ka line at 6.4 
keV. The best fit parameters for this model are given in Table|3l 



the extendend soft halo from the two regions shown in Figure^ 
in order to have the best possible signal-to-noise. 
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Fig. 6. The 0.3 - 8 keV ACIS-S spectrum of the inner 10" re- 
gion fitted with two thermal plasma components, a highly ab- 
sorbed power law and a narrow Fe Ka line at 6.4 keV. The best 
fit parameters are given in Table |3 



4. 1. The inner 10" region 

We apply to the Chandra spectrum of the inner region (10" ~ 
7.5 kpc surrounding the AGN) both models successfully tested 
in Section 3 for the XMM-Newton data. 

The thermal emission model, which was not tested in 
Paper I, when applied to the Chandra spectrum of the "nuclear" 
region, gives also a good fit (x^/dof = 65.0/75). The best fit 
parameters are given in Table |3 In this case, only two colli- 
sionally ionized plasma components are required, as shown in 
Figure|6l 

The photoionized gas model, which we discuss here, is sim- 
ilar to the one presented in Paper I. It consists of a highly ab- 
sorbed (primary) plus a less absorbed (secondary) power law 
having the same spectral index, two Fe Ka lines at 6.4 and 
6.7 keV and a single thermal plasma component. This model 
also provides a good fit (x^/dof = 62.0/75, see Figure^}. A 



column density of Nh = 6.8^2 8 ^ cm" obscuring the sec- 
ondary power law is needed, in agreement with Paper I and the 
XMM-Newton results. 

A narrow, unresolved, Fe Ka line is detected at about 3cr 
significance in the ACIS-S spectrum. Its centroid is found at 
6.34 + 0.04 keV, which is slightly lower than the expected value 
of 6.4 keV from neutral iron, probably due to uncertainties in 
the calibration. No ionized iron lines are required by the data; 
however, we derive upper limits on the fluxes of the Fe XXV 
and Fe XXVI Hnes (respectively < 0.37 x 10 and < 0.39 x 
lO"'-' erg cm"^ s"'), which are consistent with those measured 
by XMM-Newton. 

The best fit parameters for both models (summarized in 
Table ^ are consistent, within the uncertainties, with those 
found by XMM-Newton. The only exception is in the value 
of the column density obscuring the primary radiation, which 
is slightly different between the two observations (41+6 and 



Table 4. Chandra ACIS-S spectral fitting results when the ther- 
mal emission model (Model 1) and the photoionized gas model 
(Model 2) are applied to the spectrum of the inner 10". The last 
column refers to the spectrum of the extended hot gas halo in- 
stead. Fluxes are given in units of lO"'"* erg cm"^ s"', luminosi- 
ties in units of 10^^ erg s ' ; both are corrected for absorption. * 
denotes a fixed parameter. 
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servation respectively). 



4.2. Ttie liot gas tialo 

The ACIS-S spectrum of the extended hot gas halo, which 
is extracted from the region defined in Section 2.2, is source 
dominated only in the 0.5 - 0.9 keV band. Given the limited 
amount of photons (315) detected in this band, we prefer to 
use Cash statistics on the unbinned spectrum. As displayed in 
Figure |8j the spectrum is consistent with a single temperature 
plasma component with kT = 0.38 + 0.04 keV. We performed 
a Monte Carlo calculation to test the goodness of the fit and we 
found that 63% of the lO'* simulated spectra have a fit statis- 
tics less than that for the data. A fraction of approximately 
50% is expected, if the observed spectrum was produced by the 
model. If the;^'^ statistics is applied to the spectrum grouped to 
have a minimum of 20 counts per bin, fixing the model pa- 
rameters to those previously found (see Table 0}, we obtain 
do f ^15/11. 

In this analysis we adopted a more conservative approach 
in subtracting the background in comparison to Paper I. This 
resulted in a slightly smaller temperature for the extended hot 
gas halo with respect to the previous analysis, but in very good 
agreement with the XMM-Newton results. 
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Fig. 7. The 0.3-8 keV ACIS-S spectrum of the inner 10"region 
fitted with a single thermal plasma component, a highly ab- 
sorbed and a less absorbed power law with the same spectral 
index, plus a narrow Fe Ka line at 6.4 keV. The best fit param- 
eters are given in Table|3 
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Fig. 8. The 0.5-0.9 keV unbinned spectrum of the extended hot 
gas halo. A single APEC model with a temperature of 0.38 + 
0.04 keV provides a satisfactory fit. 



5. Discussion 

A broad Fe Ka emission line with cr - 0.26^q j^ keV is de- 
tected at high significance in the broad band (0.3 - 10 keV) 
EPIC-PN spectrum (Ax^ - 1 1 .4 for 3 additional free param- 
eters, corresponding to more than 99% confidence level). The 
hypothesis of a broad iron line originating from the disc is plau- 
sible, but the high column density, obscuring the AGN up to 
about 4 keV, makes the detection of the redshifted wing of a 
relativistic iron line profile more difficult. The most likely ex- 
planation for such a broad line profile is the superposition of 
a narrow, unresolved, neutral iron line at 6.4 keV, probably as- 
sociated with the AGN torus, and a blend of Fe XXV and Fe 
XXVI Unes, associated with a hot plasma, in which this tran- 
sitions are either thermally excited or induced by photoioniza- 
tion. Unfortunately, the modest statistics do not allow to sin- 
gularly resolve the three possible lines. However the fluxes and 
equivalent widths measured by Chandra and XMM-Newton are 
consistent within the two observations. 



There is no evidence for extra absorbing material obscuring 
any of the thermal plasma components. This could be due to a 
non-spherical geometry of the gas surrounding the nuclear re- 
gion, given the complex morphology of this merger, or it could 
be simply due to the low statistics. A slightly different column 
density is found to be obscuring the primary radiation between 
the two observations. Again this could be, somehow, associated 
to the complex morphology of the merger, since, for instance, 
the projected separation between the two nuclei is ~ 1", or it 
could be only an effect due to the low statistics available at high 
energies. 

According to the thermal emission model, the tempera- 
tures of the three plasma components found in Mrk 273 are 
remarkably similar to those observed in other ULIRGs (e.g. 
NGC 6240 an d Arp 22 0) and in the local starburst galaxy 
NGC 253 cPietsch et al.ll20011 The thermal emission at ~ 0.3 
keV, as clearly shown from the Chandra image, extends on 
a large scale (~ 45 kpc) and embraces the long tidal tail of 
the merger. Therefore, it must be associated with hot gas dis- 
tributed in the halo. Interestingly, thermal emission at ~ 0.7 
keV seems to be ubiquitous in the spectra of ULIRGs, prob- 
ably being associated with a nuclear or circumnuclear star- 
burst ( Fra nceschini et al1l2003b . The presence of a high tem- 
perature (> 5 keV) thermal component is less frequently ob- 
served, but not so unusual, being also detected in the XMM- 
Newton spectra of NG C 6240 (Bofler et al. 2003) and Arp 220 
jlwasawa et alJ ^2005). The ratio of the 2-10 keV unab- 
sorbed X-ray luminosity of the highest temperature ther- 
mal component to the bolometric luminosity in Mrk 273 is 
L2-iQk,v /Lfir ^1x10-^ {Lfir = m - lOOOyum ~ Lh„i), a 
factor of 20 smaller than in NGC 6240 and a factor of 5 larger 
than in Arp 220. 

On the other hand, the second model investigated suggests 
that a secondary hard X-ray component originates from the re- 
flection of the primary continuum on some photoionized cir- 
cumnuclear materia l. This is often observed in Compt on-thick 
Seyfert 2 galaxies teambr unaeta L.2Q01: KinkhabwalaetalJ 
l2002t iDone et alJl2003t iBianchi et alJl2005h . which, in some 
cases, also show evidence for starburst activity (e.g. Circinus, 
Ruiz et al. 2001; NGC 1068, Smith & Wilson 2003). As was 
observed in Paper I, two bright clumps of soft X-ray emission 
coincide with [OIII] emission. This piece of evidence seems to 
support the hypothesis of a photoionized gas surrounding the 
nucleus, since in several Compton-thick Seyfert galaxies, the 
extended soft X-ray emission expected to be produced by re- 
flection from that photoionized gas has been observed to coin- 
cide with the [OIII] emission and the ionization cone (e.g. Sako 
et al. 2000; Young et al. 2001). 

The absorption corrected hard X-ray luminosity of 
Mrk 273 is L2-iokeV ~ 7 x lO'*^ erg s"', which corresponds to 
~2xl0"^ of the far-IR luminosity (L/r/R ~ lO'^L©, e.g. Genzel 
et al. 1998). Such a ratio approaches values of pure starbursts 
and is orders of mag nitude smaUer than in Compton-thin AGNs 
(Risaliti et al. 2000). According to the thermal emission model, 
the absorption coiTected 0.3 - 10 keV luminosity due to the 
three thermal components would be approximately 0.7 x 10"*^ 
erg s"'. Comparable soft X-ray luminosities have been mea- 
sured for NGC 6240 (iNetzer et aL.2005.> and other ULIRGs in 
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which the starburst contribution t o the b olometric luminosity 
is dominating ( Franceschini et al] l2003l) . From the photoion- 
ized gas model the resulting thermal contribution to the X-ray 
luminosity is about a factor of 3 smaller For Mrk 273 the un- 
obscured AGN X-ray luminosity is a more modest fraction of 
the bolometric luminosity than in typical AGN. This suggests 
additional strong star formation which will naturally explain 
the thermal X-ray components invoked in one of our two sce- 
narios. The presence of strong star format ion agrees with the 
detection from mid-infrared spectroscopy jGenzel et alJll99^ 
of both very powerful star formation and an AGN in Mrk 273. 

It must be mentioned that, besides the two models dis- 
cussed above, other interpretations are possible. In fact X- 
ray emission fro m luminous radio supernovae like SN1995J 
( IFox et a l. 2000), or the integrated contribution from low- and 
high-mass X-ray bin aries, as observed from t he spectra of lo- 
cal starburst galaxies jPersic & Rephaell2002l) . may both give 
rise to a very hot (kT ~ 5 keV) t hermal component, or a F ~ 2 
power law. Following lKennicut3 (>199&) . the star formation rate 
estimated assuming that the entire Lfir is due to starburst 
would be ~ 200 M© yr Therefore, the inferred supernova 
rate is ~ 2 yr which i s in agreement with the value estimated 
from radio observations lBondi et al l ll2005l) . On the other hand, 
given a 200 M© yr"' star formation rate, the e stimated emis- 
sion from X-ray binaries, following iFranceschini et aL C2003 ). 
would be L2-iokeV ~ 2x10"*' erg which is very close to the 
observed value (~ 3 x lO"*' erg s"'). Finally, it is worth noting 
that X-ray spectra from supernovae have shown evidence for 
strong Fe K line at 6.7 keV (e.g. SN1986J, Houck et al. 1998), 
therefore could also account for this feature, while emission 
from X-ray binaries can not. 

6. Conclusions 

We analyzed the X-ray spectral properties of Mrk 273 combin- 
ing the high throughput of XMM-Newton with the high spatial 
resolution of Chandra. 

From the XMM-Newton spectrum we found a broad Fe Ko- 
emission line with a high statistical significance (>99% c.l.). 
We suggested the superposition of multiple unresolved iron line 
features: one from neutral iron at 6.4 keV {EW ~ 170 eV) and 
one from a blend of Fe XXV at 6.7 keV (EW ~ 120 eV) and 
Fe XXVI at 6.97 keV (EW < 85 eV). 

We tested two different models capable of reproducing the 
possible ionized iron lines detected: one through thermal emis- 
sion, therefore involving the presence of a hot gas associated 
with a starburst, and the other one through reflection of the 
primary continuum by some optically thin photoionized gas. 
Unfortunately, given the available statistics, we can not rule 
out any of the two possible explanations. Some indications that 
a reflection componet is indeed present in the nuclear region 
of Mrk 273 may come fro m polarization measures, which ar e 
already available in the IR jSiebenmorgen & Efstathioul200ll) . 
but still not in the X-rays. High resolution spectroscopy could 
be helpful in order to detect some of the forbidden lines or the 
Radiative Recombination Continua (RRC), which are consid- 
ered as typical signatures of photoionization ('P orquet & DubaJ 
l2000^.Liedahn999:.Liedahl & Paerels.1996.) . 
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Fig. A.l. The combined MOS 1 and MOS2 0.3 - 10 keV data 
fitted with a power law modified by Galactic absorption. 

Appendix A: Thie unabsorbed Seyfert 2 
Mrl< 273x 

Over recent years, much attention has been directed to the study 
of unabsorbed Seyfert 2 galaxies (e.g. Ptak et al. 1996; Panessa 
& Bassani 2002; Gallo et al. 2005). These galaxies are of in- 
terest for two reason: they may account for as much as 30% of 
the Seyfert 2 population (Panessa & Bassani 2002) . and they 
are a challenge to explain with the existing Seyfert unification 
model. 

Mrk 273x (z - 0.458) is one of the better known unab- 
sorbed Seyfert 2 galaxies, perhaps because it falls in the field- 
of-view of most observations of Mrk 273. In this appendix, we 
would like to discuss the properties of Mrk 273x as seen during 
the XMM-Newton observation of Mrk 273. 

In the EPIC-PN image, Mrk 273x falls directly in a chip 
gap; hence we used the combined MOS 1/2 data for spectral 
analysis. The combined spectrum is source dominated in the 
0.3 - 10 keV band (0.44 - 14.6 keV in the rest-frame). The 
total number of source plus background counts is 526, whereas 
the total number of background counts is 65 (scaled to source 
cell size). 

The best-fit continuum model of Paper I to the Chandra 
data of Mrk 273x was a power law plus a weak thermal compo- 
nent. Significantly detected above the continuum was an emis- 
sion line at 1.82 keV, which the authors suggested could be due 
to Si XVI. The XMM-Newton data are consistent with only a 
power law (F - 1.49 ±0.12) modified by Galactic absorption 
(X^/dof = 60.2/82; FigureOl. 

The 90% confidence upper limit on the intrinsic absorption 
is < 4.5 X lO^'' cm"^, at least a factor of three less than what 
was required in the Chandra model. 

Multi-component continuum models (e.g. blackbody plus 
power law, double power law, or a broken power law) also fit 
the data well, but they are not significant improvements over 
the single power law. Neither the addition of a ~ 1.82 keV 
emission feature is required by the data (Ax^ = 1 for 2 ad- 
ditional free parameters; EW < 33 eV). Positive residuals at 
about 4.3 keV can be modelled with a narrow Gaussian profile 
with a rest-frame energy of E - 6.29;^° 34 keV. The addition 
of the Gaussian profile is an improvement to the residuals, but 
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it is not statistically significant (Ax^ - 3 for 3 additional free 
parameters). 

Based on the single power law fit, the unabsorbed 0.5 - 
10 keV flux is (1.43+JJ-22) x 10"'^ erg s"' cuT^ (1 .18 x lO"'^ erg 
s"' cm"^ in the 2-10 keV band). The intrinsic 0.3 - 10 keV 
unabsorbed luminosity is 9.2 x lO"*"* erg s"'. The broad band 
flux measured during this XMM-Newton observation is com- 
parable to the flux observed during the Chandra observation 
two years earlier. The ~ 22 ks light curve of Mrk 273x shows 
no variability and is perfectly fitted by a constant. 
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